Stability analysis of singular patterns in the 1 - D Gray {

In this work, we analyze the linear stability of singular homoclinic stationary solutions and spatially-periodic stationary solutions in the one-dimensional Gray-Scott model. This stability analysis has several implications for understanding the recently discovered phenomena of self-replicating pulses. For each solution constructed in 5], we analytically nd a large open region in the space of the two scaled parameters in which it is stable. Speciically, for each value of the scaled inhibitor feed rate, there exists an interval, whose length and location depend on the solution type, of values of the activator (autocatalyst) decay rate for which the solution is stable. The upper boundary of each interval corresponds to a subcritical Hopf bifurcation point, and the lower boundary is explicitly determined by nding the parameter value where the solution`disappears,' i.e., below which it no longer exists as a solution of the steady state system. Explicit asymptotic formulae show that the one-pulse homoclinic solution gains stability rst as the second parameter is decreased, and then successively, the spatially periodic solutions (with decreasing period) become stable. Moreover, the stability intervals for diierent solutions overlap. Explicit determination of these stability intervals plays a central role in understanding pulse self-replication. Numerical simulations connrm that the spatially periodic stationary solutions are attractors in the pulse-splitting regime; and, moreover, whenever, for a given solution, the value of the activator decay rate was taken to lie in the regime below that solution's stability interval, initial data close to that solution was observed to evolve toward a diierent spatially periodic stationary solution, one whose stability interval included the parameter value. The main analytical technique used is that of matched asymptotic expansions.